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Abstract 

Electromagnetic properties of hadrons can be computed by lattice simu- 
lations of QCD in background fields. We demonstrate new techniques for 
^ • the investigation of charged hadron properties in electric fields. Our current 

O i calculations employ large electric fields, motivating us to analyze chiral dy- 

namics in strong QED backgrounds, and subsequently uncover surprising non- 
perturbative effects present at finite volume. 

> 

(N : 1 Introduction and Motivation 

Electromagnetic properties paint an intuitive picture of the distribution of charge and 
Q ■ magnetism within a hadron. Electromagnetic polarizabilities, for example, encode the 

On . ability of a hadron to deform in response to an applied field. This stiffness described by 

^ ' polarizaiblities is a fundamental property of the hadron, with a transparent physical 

meaning. Chiral dynamics, moreover, places tight constraints on the polarizabilities, 
^ ' which, in some cases, have not compared well with experiment. 

■ Electromagnetic properties of hadrons can be computed by lattice simulations of 

QCD in background fields. For electromagnetic polarizabilities, the background field 
method is currently the only option because computations using four-point functions 
are not feasible. Our lattice QCD computations in background fields and related 
theoretical developments are summarized here. We demonstrate new techniques for 
the investigation of charged hadron properties in background electric fields. Our cur- 
rent calculations employ large electric fields, motivating us to analyze chiral dynamics 
in strong QED backgrounds, and subsequently uncover surprising non-perturbative 
effects present at finite volume. 



^Talk given by B. C. Tiburzi at the Tenth Workshop on Non-Perturbative QCD at I'lnstitut 
d'Astrophysique de Paris, France, 8-12 June 2009. 
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2 Background Field Method 



There are two basic steps involved in background field computations of hadron prop- 
erties. First one computes hadron correlation functions using lattice QCD in classical 
external fields, for various values of the external field strength. One then compares the 
calculated correlation functions with their behavior predicted by the relevant single- 
particle effective action. The parameters of the single-particle effective action can be 
deduced by considering the field-strength dependence of the correlation functions. 

The simplest application of the method is for a neutral scalar particle. We use 
the neutral pion to exemplify this. Imagine correlation functions for the neutral pion 
have been calculated in a background electric field specified by the gauge potential 
Af^{x) = —Sx4^6^^3. As a composite scalar, the most general effective action for the 
neutral pion in an external electric field has the form 

£(p = 0) = ^7r° -8484 + mlo +171^0 a eS'^ 7r°. (1) 

Here we work in Euclidean space, and have projected the Lagrangian onto vanishing 
spatial momentum, which is possible due to our implementation of the electric field. 
Terms involving higher powers of the electric field are possible, but we have kept 
only the lowest-order contribution as we imagine the fields are sufficiently weak. The 
coefficient of the second order term, a^;, is the electric polarizability. 

Using Eq. ([1]), one can deduce the behavior of the neutral pion two-point correla- 
tion function, namely 

G(r) =^(0|7rO(f,r)7rO(0,0)|0) = Zexp(-Er), (2) 

where the neutral pion energy, E, is given by E = m^o + ^aE^"^, and the sign arises 
from our treatment in Euclidean space. Thus by determining the neutral pion energy 
as a function of the applied electric field, one can deduce the electric polarizability. 

Recently we have carried out this procedure for the neutral pion (additionally for 
the neutral kaon) [Tj. We used an ensemble of anisotropic (2 + l)-fiavor dynamical 
clover lattices generated by the Hadron Spectrum collaboration |2l [3]. Specifically 
our ensemble consists of 200 configurations of size 20^ x 128 with light quark masses 
giving rise to a pion mass of 390 MeV. The anisotropy factor has been tuned to 
O's/O't = 3.5. We computed neutral pion correlation functions for various values of 
the electric field, and then fit to the expected form in Eq. ([2]). Due to our choice of 
lattices, we use a two-state fit in order to stabilize the extraction of the ground-state 
energy. In Figure [1], we show one such fit for the lowest value of the electric field 
(denoted by n = 1). Additionally we show the energies extracted as a function of 
the external field strength. For sufficiently small fields, the behavior of the energies 
should be quadratic with respect to the field. The fit shown (denoted II) includes 
both quadratic and quartic terms. 
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Figure 1: On the left: fit to the neutral pion correlation function in the electric 
field with strength corresponding to n = 1. Plotted is the effective mass Me//(t) = 
— \ogG(t + 1)/G(t). On the right: a fit (denoted II) to the electric-field dependence 
of the neutral pion energy. 



3 Field Strengths and Strong Fields 

Having investigated the neutral pion energy as a function of the external electric field 
strength, we notice that the values of the external field employed are not necessarily 
perturbatively small. We are not able, however, to simulate at arbitrary values of the 
field. The implementation of a uniform field coupled to matter fields is restricted by 
quantization conditions relevant for a torus [1] . 

For our background electric field, the field strength £ must be quantized in the 
form 

e = (3) 

where qa = — l/3e is the electric charge of the down quark, L is the length of the 
lattice in the Xa-direction, T is the length in the a;4-direction, and n is an integer. To 
implement a uniform field on a torus, one must use link variables. The standard links, 
U^\x) = exp[iqAfj_{x)] = exp^—igSx^S ^^3) , are not enough to ensure a uniform field 
on a discrete torus. Plaquettes that wrap around the time-boundary see a dramatic 
spike in the electric field. This effect can be removed by quantizing the fields, and 
adding additional transverse gauge links U^{x) = exp{iq£Tx36fj,^4^6x4^,T-i) \5\- The 
net result is to modify the color links U^{x) in the following way 

U,{x)-^U,{x)U'^^{x)Ui{x). (4) 

Notice this post-multiplication leaves the sea quarks electrically neutral. Charging 
the sea would require at least an order of magnitude greater computing resources. 

To investigate the necessity of quantized field strengths, we plot in Figure [2] the 
neutral pion correlation function calculated in two different classical fields. Both 
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7r°: Eq. (4), tsrc=52 
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Figure 2: Comparison of quantized and non-quantized field strengths. In the plot, 
correlation functions have been translated forward by twelve units in time. 



fields are implemented using Eq. (jlj), but we use the values n = e = 2.71828 . . ., and 
n = 3. While numerically the fields are comparable in magnitude, the non- integer 
value results in a single plaquette at the far corner of the lattice with a dramatic 
non-uniformity. We can see the effect of the non- uniformity by letting the neutral 
pion propagate around the time boundary. For this exploration, we use configurations 
with 64 sites in the time direction, and place the pion source at t = 52 (details of this 
study were presented in [6]). After translating the correlation functions forward by 
twelve units, the onset of a plateau in the effective mass at around t = 12 corresponds 
to the edge of the lattice in the time direction. A plateau is only seen for the n = 3 
value of the field. 

While current lattice sizes give rise to rather large values of the external fields, we 
are still in the regime of a so-called strong field power counting in chiral perturbation 
theory. On the 20^ x 128 anisotropic lattices, we have {e£ /mlY = 0.18n^, where n is 
the integer in the field strength quantization condition. For such values of the field, 
we can perform chiral perturbation theory calculations to determine the response 
of the pion cloud to the strong field [7]. Such calculations treat chiral corrections 
perturbatively m^/A^ ^ 1, but re-sum charge couplings to all orders, e£/ml ~ 1. 

As an example of a strong field calculation in chiral perturbation theory, one can 
derive the field-strength dependence of the neutral pion energy in strong electric fields. 
The shift, 6E, of the neutral pion energy has the form 



ml 



E{p = 0) — m^r 




m^ 

where the explicit expression for I{x) arising at one-loop order has the form 



(5) 



T(x) = X ( 1 - log + 2 logF (^i^ ) - log27r. (6) 
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Figure 3: Plot of the electric field dependence of the neutral pion energy. The dashed 
curve shows the non-perturbative result at one-loop order in the strong field chiral 
expansion. LO and NLO represent the leading and next-to-leading order perturbative 
expansions of the energy in powers of the electric field, where LO oc and NLO 
includes terms oc S"^. Successive perturbative approximations are also depicted. 



In Figure [3l we plot the neutral pion energy as a function of the electric field, ^ = 
eS/V6ml. The non-perturbative result of Eq. ([5]) is contrasted with the first two 
non- vanishing perturbative approximations. In the range of ^ shown, the perturbative 
expansion has broken down: higher order perturbative approximations result in worse 
agreement with the non-perturbative result. Consequently the LO approximation 
does best. This behavior is characteristic of asymptotic expansions. 

Finally notice that on the same lattices, magnetic fields are further from the per- 
turbative regime, i.e. we have {eB/mlY = / L"^ {e S / , due to the quantization 
condition, with T'^/L'^ = 3.3 for our anisotropic lattices. 



4 Charged Pion 

The background field method can also be applied to charged particles, as was first 
suggested in [8]. The philosophy is the same: compute the hadron correlation function 
in various values of the external field strength, and then match onto the behavior 
expected from the single-particle effective action. 

Let us use the charged pion as an example. Imagine we have computed the charged 
pion correlation function for various values of an external electric field. What should 
we expect for the behavior of these correlation functions? For the charged pion, the 
single-particle effective action has the form 

C{p = 0) = TT^ [-94^4 + S^xl + E\S)\ TT. (7) 

The explicit time-dependent term arises from minimal coupling of the electric field 
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to the charged pion. The non-minimal couphngs have been subsumed into the rest 
energy, E{S), which has a small-field expansion in powers of the electric field. The 
coefficient of the first field-strength dependent term in this expansion is proportional 
to the polarizability. 
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Figure 4: Charged pion correlation functions in external electric fields. 



Using Eq. ([7]), we deduce the form of the charged pion correlation function, namely 
^(r) =5^(0 1 7r+(f,r)7r- (0,0)1 0) = ZVe^ (V2S^) , (8) 

X 

where T>i^{z) is a parabolic cylinder function (for which a convenient integral rep- 
resentation can be derived using Schwinger's proper-time trick [SIE])- In Figure HI 
we show the behavior of the charged pion correlation function for four values of the 
external field strength, corresponding to the integers n = 1-4. Shown along with the 
correlation functions are fits using the relativistic scalar propagator, Eq. ([8]). As with 
the neutral case, we have performed a two-state fit. The anisotropy is of considerable 
help yielding 3.5 times the number of time slices over which to fit the correlation 
function. While the pion charge leads to correlation functions with dramatically dif- 
ferent time-dependence than a simple exponential falloff, the behavior seen is inline 
with expectations for a relativistic charged scalar in a constant field. 
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5 Finite Volume Artifacts 



While the lattice volume restricts the strength of external fields, there are additionally 
finite size effects of a dynamical nature. Physically one expects rather sizable modifi- 
cation to the electromagnetic polarizabilities in finite volume because these quantities 
physically have the dimension of volume. Furthermore a hadron's deformation arises 
from its charged pion cloud, and lattice volumes are not considerably greater than 
the pion Compton wavelength, Ac = I/jtItt ~ 3/2 fm. It is naive, however, to merely 
ascertain how much of the induced dipole moment owes to finite volume [ini flT] . 
As multipole polarizabilities arise from 5*0(3) invariant physics, we must instead 
deal with the subtleties of electromagnetic fields on a torus, which only maintains a 
discrete cubic subgroup of rotations. 

A related subtlety in the physics on a Euclidean torus is the lack of boost invari- 
ance. Matrix elements in the rest frame are related to those in a boosted frame only 
up to a finite volume modification. Thus the rest frame matrix element for the pion 
charge 

(7r(0)|J4|7r(0)) = 2p4Q (9) 
does not generate the current J = Qv when boosted to an arbitrary frame, 

{nip)\JHp)) = 2p[Q - Ag(L)]. (10) 

This fact is consistent with Ward-Takahashi identities, and is analogous to charge 
screening at finite temperature [12]. While Ward identities are violated, they are 
not applicable in finite volume because a limiting procedure does not exist to take 
derivatives with respect to the quantized momentum modes, p = 27in/L. One can use 
coordinate space methods in finite volume to derive finite volume corrections using 
chiral perturbation theory \X3\. Specifically one can derive the single-particle action 
in finite volume; which, for the charged pion, has the form 



AC{p = 0) = vrt \AQ\L)S^xI + Aa{L)S' 



vr, ill] 



where only terms up to 0{£'^) have been retained. 

To our surprise, the perturbative expansion of finite volume corrections is not 
sufficient to account for all of the systematics present in the lattice data. Shown in 
Figure [5] is an unexpected systematic effect: the extracted rest energy of the charged 
pion shows sensitivity to the fit window chosen. Finite size effects can give rise to 
such oscillatory time dependence as we now argue. 

Using our gauge potential for the electric field, there are other gauge invariant 
quantities on a torus besides the field strength £. These quantities are Wilson loops 
formed using the gauge holonomy, B, given by 

Q = dxsA^ix) = -SLx^. (12) 
Jo 
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Figure 5: Extracted charged-pion rest energies as a function of the fit window for the 
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3 field strength. The fit window spans [tmin, tmax], with tmax = 30 fixed. 



Thus the Wilson loops, W{Q) = e*®, will be time dependent, and such oscillatory 
behavior cannot be captured in perturbation theory in the strength of the field. One 
needs a finite volume calculation that treats 0(x4) ~ 1. Results of these chiral 
perturbation theory calculations are rather cumbersome to work with, and we can 
make considerable progress by appealing to symmetries. 

Consider the neutral pion effective action in finite volume. The operators appear- 
ing are constrained by parity, charge conjugation, gauge invariance, and discrete time 
translation. Writing down all terms consistent with these symmetries we find. 



C{p = 0) = -n' 



n=l 



TT 



(13) 



where all parameters implicitly depend on the lattice size L. The relevant Wilson 
loops, >V(x), are Hermitian, and given by >V(x) = |[iy(x) + iy^(a;)]. Physically 
the tower of Wilson loop couplings arises from charged virtual pions propagating 
around the lattice n times. Given that the probability for such propagation scales 
exponentially, Vn oc exp(— nmTr-L), it is reasonable to consider only the first term in 
the tower of Wilson loop operators. The resulting single-pion effective action 



-TT 



TT 



(14) 



leads to a pion propagator that can be expressed in terms of Mathieu functions. It 
seems likely that such oscillatory functions can explain the variation of the extracted 
energies with respect to the fit window. This work, however, is still in progress. 
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6 Summary 



Above we have reviewed lattice QCD computations in background fields. We have 
stressed the two basic features of such calculations; namely, the measurement of 
hadronic correlation functions in external fields, and the matching of the expected 
behavior for these correlation functions obtained from single-hadron effective actions. 
For ease we only considered pions in background electric fields. 

The neutral pion in infinite volume presents the simplest application of the back- 
ground field method: a field-strength dependent shift of the pion energy is the only 
modification to the neutral pion correlation function. Measuring this shift allows 
one to deduce the electric polarizability. The charged pion electric polarizability can 
also be deduced by using the background field method. For this case, however, the 
calculated correlation function has non-standard behavior as a function of time. Nev- 
ertheless, this behavior of the correlation function is precisely that predicted for a 
relativistic charged particle in a constant electric field. 

Our computations are currently limited by systematic errors. As field quantization 
conditions must be met, the relatively small lattice volume employed leads to large 
values of the electric field. As one leaves the regime in which the field-strength 
dependence is perturbativc, strong field chiral perturbation theory can be applied. In 
fact, chiral perturbation theory gives parameter- free predictions for hadron energies 
at next-to-leading order in the chiral expansion as a function of eS/m^. Small lattice 
volumes lead to another malady: finite size effects. Single pion effective actions are 
more complicated in finite volume. Quantum fiuctuations allow virtual pions to wrap 
around the lattice and encounter background fields that differ topologically. Couplings 
of Wilson loops to pions must be included in the single pion effective action, and can 
potentially explain the oscillatory behavior seen in extracted energies as a function 
of the fit window. 

Lastly, our computations have been carried out at a single value for the pion mass, 
and with vanishing electric charges for the sea quarks. Further studies are required 
to extract physical results. Nevertheless, we have demonstrated new techniques for 
background field simulations which warrant such further studies. We additionally 
intend to study spin-1/2 baryons using generalizations of the methods presented here. 
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Discussion 

J. Rafelski (University of Arizona): Throughout you have referred to the neutral 
pion. Is this correct? What about isospin? 

Tiburzi: In background electromagnetic fields, isospin symmetry can no longer be 
used to simplify the calculation of neutral pion correlation functions. We have only 
calculated the connected part of the correlation function, and so "neutral pion" should 
be thought of as a convenient name rather than the physical propagating state. The 
calculation of the disconnected part is left for future investigation. 
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